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Abstract 

Conformally flat spherically symmetric cosmological models rep¬ 
resenting a charged perfect fluid as well as a bulk viscous fluid dis¬ 
tribution have been obtained. The cosmological constant A is found 
positive and is a decreasing function of time which is consistent with 
the recent supernovae observations. The physical and geometrical 
properties of the models are discussed. 
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1 Introduction 


A considerable interest has been shown to the study of physical properties 
of spacetimes which are conformal to certain well known gravitational fields. 
The general theory of relativity is believed by a number of unknown func¬ 
tions - the ten components of g l K Hence there is a little hope of finding 
physically interesting results without making reduction in their number. In 
conformally flat spacetime the number of unknown functions is reduced to 
one. The conformally flat metrices are of particular interest in view of their 
degeneracy in the contex of Petrov classification. A number of conformally 
flat physically significant spacetimes are known like Schwarzschild interior 
solution and Lemaitre cosmological universe. 

At the present state of evolution, the universe is spherically symmetric 
and the matter distribution in it is isotropic and homogeneous. Buchdahl 
[1] has obtained the conformal flatness of the Schwarzschild interior solution. 
Singh and Roy [2] have discussed the possibility of existence of electromag¬ 
netic fields conformal to some empty spacetime. Singh and Abdussattar [3] 
have obtained a non-static generalization of Schwarzschild interior solution 
which is conformal to flat spacetime. Roy and Bali [4] have obtained the 
solution of Einstein’s field equations representing non-static spherically sym¬ 
metric perfect fluid distribution which is conformally flat. Pandey and Tiwari 
[5] have discussed conformally flat spherically symmetric charged perfect fluid 
distribution. Reddy [6] and Rao and Reddy [7] discussed static conformally 
flat solutions in the Brans-Dicke and Nordtvedt-Barker scalar-tensor theo¬ 
ries. Shanthi [8] has shown that the most general conformally flat static 
vacuum solution in the Nordtvedt-Barker scalar-tensor theory is simply the 
empty flat spacetime of general relativity. There has been a recent literature 
(Mclfo and Rago [9], Mannheim [10], Yadav and Prasad [11], Endean [12,13], 
Obukhov et al. [14], Mark and Harko [15]) which shows a significant interest 
in the study of conformally flat spacetime. 

Most cosmological models assume that the matter in the universe can 
be described by ‘dust’(a pressureless distribution) or at best a perfect fluid. 
Nevertheless, there is good reason to believe that - at least at the early stages 
of the universe - viscous effects do play a role [17]-[19]. For example, the ex¬ 
istence of the bulk viscosity is equivalent to slow process of restoring equilib¬ 
rium states [20]. The observed physical phenomena such as the large entropy 
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per baryon and remarkable degree of isotropy of the cosmic microwave back¬ 
ground radiation suggest analysis of dissipative effects in cosmology. Bulk 
viscosity is associated with the GUT phase transition and string creation. 
Thus, we should consider the presence of a material distribution other than 
a perfect fluid to have realistic cosmological models (see Gr0n [32]) for a re¬ 
view on cosmological models with bulk viscosity). The effect of bulk viscosity 
on the cosmological evolution has been investigated by a number of authors 
in the framework of general theory of relativity [21]-[31]. 

The purpose of this paper is to apply conformally flat spherically sym¬ 
metric line element to charged perfect fluid and to bulk viscous fluid models 
in cosmology. This paper is organized as follows. The held equations are 
presented in Section 2. Section 3 includes the solution of the held equations 
in presence of charged perfect fluid distribution. Section 3.1 contains some 
physical properties of the model. Finally in Section 4 the bulk viscous models 
are considered. 


2 Field Equations 

We consider the conformal metric in spherical polar coordinates 


ds 2 = e x (dr 2 + r 2 dd 2 + r 2 sin 2 6d(j) 2 — dt 2 ), (1) 

where A is a function of r and t. We number the coordinates as x l = r, 
x 2 = 9, x 3 = (J) and x A = t. 

The energy momentum tensor for distribution of a charged perfect fluid has 
the form 

Tij = (e + p)viVj + pgij + E ih (2) 

where Eij is the electromagnetic held given by 


Eij = — 

v Att 


FaiF bj g ab - gijF ab F 1 


ab 


( 3 ) 


Here e and p are the energy density and isotropic pressure respectively and 
v l is the how vector satisfying the relation 


9%jv'v? = - 1 . 


( 4 ) 
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The electromagnetic field tensor F tJ satisfies Maxwell’s equations 


F l j = 4npv\ (5) 

F[ij-,k] = 0 , ( 6 ) 

p being the current density. Here and henceforth a comma and a semicolon 
denote ordinary and covariant differentiation respectively. The Einstein field 
equations 


Fjj g 9ijF T b (j ij SirTij, 

for the line element ( 1 ) has been set up as 

r/ N o 3A 2 2Ai , \\ 

8 yr[(e + p)v{ + pe x ] = — + — - A 44 - — 


+ e _A (F 14 ) 2 + Ae, 
A? , A, . A 2 


87 t pe — An H—■—I- A44 --— e (F14) + Ae , 

4 r 4 

2 „„ai_3A1_ a _A|_2Ai 


8n[(e + p)vl -pe x ] = 


4 4 

-e- A (Fi 4 ) 2 -Ae A , 

87r(e + p)v 1 U 4 = —^ - A 14 . 


r 


(7) 


( 8 ) 

(9) 


( 10 ) 

( 11 ) 


Equation (|D gives 



( 12 ) 


3 Solutions of the field equations 

From eqs. (||) and (|9|) we have 

8' 7r [( e + p) v l] ~ 2e a (Ai 4 ) 2 = -A -|-An. (13) 

2 r 

Also eqs. @ and (|T 0 |) readily give 

8yr[(e + p)u 2 ] + 2e _A (Fi 4 ) 2 = ^ - — - A 44 . (14) 

2 r 
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Combining eqs. (]T3|), ([13]) and we obtain 


87r[(e + p)e x ] + 4e X (F U ) 2 


A|_ A|_ 2Ai 
2 2 r 


— A44 + An. 

Equations @ and flT5[) together reduce to 

8nee x + 3e _A (iq 4 ) 2 = ^ f A 2 - A 2 - 


In comoving coordinate system v\ = 0, then eq. (|T3|) reduces to 

_ = M + h. _ hi. 

V ’ 4 2r 2 

From eq.(|Tl|) we obtain 

2A14 — A1A4 = 0. 

The general solution of (18) is obtained as 

e x = [a(r) + /3(«)]' 2 , 


(15) 

(16) 

(17) 

(18) 
(19) 


where a and (3 are functions of r and t respectively. 

Hence the geometry of the spacetime (1) reduces to the form 

ds 2 — 1 - (dr 2 + r 2 d6 2 + r 2 sm 2 0d0 2 — dt 2 '] , ( 20 ) 

(a + /?)- v 7 ' 

which is the model for a distribution of charged perfect fluid with the flow 
vector in t-direction. The pressure and density for the model (20) are given 
by 

( 3 \ 

2^44 - an-— J , (21) 

87 re + A = 3(/3 2 — a 2 ) + 3(cr + /?) H—-—^ . (22) 

Let us assume that the fluid obeys an equation of state of the form 


V = 7 e > 


(23) 
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where 7(0 < 7 < 1) is constant. Using eq. (|2B| ) in eqs. ( |2T| ) and eq. (^2[), we 
get 

£ = sl“l + +7) (Al + “ n) ’ (M) 

A = —j- ——^■(o: + /?) (/?44 + <Ui) — 3(ay — (3%) 

(1 + 7) 

+ (a + {3)((3 44 —2an). (25) 

If we put A = 0 in our solution, we recover the solution obtained by Pandey 
and TiwariQ. 

Particular cases: 

Case (i): If we consider (3{t ) = |r; a, a > 0, eqs. (^4j) and (|25| ) reduce to 


e = kiat 4 + kiat 6 , 

(26) 

A = + a(^’2 + 12a)f -6 , 

(27) 

where h - 4 ^ a +y) ,k : 2 - (1+7 ]. 

Case (ii): If we consider j3{t ) = f; a, a > 0, eqs. (|| 

) and (25) reduce to 

e = ak 3 t~ 3 + ak 3 t ~, 

(28) 

A = ak 3 t~ 3 + a(3a + Aq)t^ 4 , 

(29) 


where k 3 = 1?r( ,° + ^ , fc 4 = jyyyy- It is observed from eqs. (^7j) and (|29|) that 
the cosmological constant A, in both cases, is a decaying function of time and 
it approaches a small value as time progresses (i.e., the present epoch), which 
explains the small value of A at present. Additionally, A also comes out posi¬ 
tive which is consistent with the recent supernovae observations (Perlmutter 
et ah [33], Riess et ah [34], Garnavich et ah [35], Schmidt et ah [36]). 

3.1 Physical properties of the model 

The non-vanishing component of the flow vector, v 4 is given by 

1 

V4 (a + (3) 

The reality conditions (e + p) >0 and (e + 3 p) > 0 lead to 

aq 

P44 + ctii H-> 0, 


(30) 

(31) 
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and 


(32) 


(a) - 0l) + (a + P) (^44 - > 0 . 

Using eq. (|T9|) in eq. (|I7|) gives 

(v - an)' 


= 


(a + f3 ) 


From eqs. (® and (|33|) the current density p is given by 


P = 


{a + (3) 3 d 

rp2 


(v-«n) : 


(a + (3 ) 


(33) 


(34) 


The non-vanishing component of the acceleration vector 

v\ = v i . j v 3 , (35) 


is given by 


v i = 


Of i 

(of + /3) 


(36) 


Thus the acceleration is always directed in radial direction and the fluid flow 
in t —direction is uniform. If aq < 0 , acceleration is positive and if oq > 0 , 
there will be deceleration. 

The expression for kinematical parameter expansion 6 is given by 


6 = 3(3 4 , (37) 

All components of rotation Wij and shear ay,- tensors are found to be zero. 
We observe that the expansion is time-dependent only. Hence the model (20) 
representing a distribution of charged perfect fluid is expanding with time 
but non-rotating and non-shearing. 


4 Bulk viscous models 

In this section bulk viscous models of the universe are discussed. Weinberg 
[16] has suggested that in order to consider the effect of bulk viscosity, the 
perfect fluid pressure should be replaced by effective pressure p by 

p = p-£d, (38) 


7 










where p represent equilibrium pressure, £ is the coefficient of bulk viscosity 
and 9 is the expansion scalar. Here £ is, in general, a function of time. 
Therefore, from eq. (EH), we obtain 

87 r(p - £0) - A = 3(oq - £?|) 

( 3 \ 

2/?44 — Qqi-J . (39) 

Thus, given £(£) we can solve the cosmological parameters. In most of the in¬ 
vestigations involving bulk viscosity is assumed to be a simple power function 
of the energy density [21]-[23]. 


m = £oe n , (40) 

where £0 and n are constants. If n = 1, eq. (|40j) may correspond to a 
radiative fluid. However, more realistic models [24] are based on lying in the 
regime 0 < n < 4. 


4.1 Model I : (£ = £ 0 ) 

In this case we assume n — 0 in eq. ( |TU[ ) which gives £ = £0 = constant. By 
the use of eqs. (|3|) and fl37P in eqs. (^2]) and (|39|) , we obtain 


47 t (1 + 7)6 — 127t£ 0 /?4 + (cr + f3){a.n + ^ 44 ) 


(41) 


(1 + y)A — 3(1 + 7)(/?| — a'l) 


+ (o; + (3) 


Cti 

(1 + 37 ) 0(11 + 3(1 + 7) -2/^44 

r 

- 247t£o/3 4 


(42) 


4.2 Model 11 : (£ = £oe) 


In this case we assume n — 1 in eq. 
eqs. (|23[) and (|37|) in eqs. (|22|) and ( 


which gives £ 
|), we obtain 


£ 0 e. By the use of 


{a + (3) (at 11 + ^ 44 ) 

(1 + 7 — 3£ q /?4) 


(1 + 7)A — 3(1 + 7)(/3| — aq) 
8 


( 43 ) 









+ (ct + /?) 


(1 + 37)011 — 2/?44 + 3(! + 7) 
6£o/^ 44 (tt + /3)(o ; ll + ^ 44 ) 


Oi 

r J 


These eqs. (|4T|) - 
of /3(f). 


(44) 

(1 + 7 — 3^ o /34) 

will supply different viable models for suitable choices 


5 Conclusions 

We have obtained conformally flat spherically symmetric cosmological mod¬ 
els in the presence of a charged perfect fluid where the acceleration is always 
directed in radial direction and the fluid flow in f-direction is uniform. We 
have also discussed two particular cases. In both cases we observe that the 
energy conditions hold good and the cosmological constant is found posi¬ 
tive and is decreasing function of time which is consistent with the present 
observations. The model is expanding with time but non-rotating and non¬ 
shearing. 

Assuming £(f) = £ 0 e n , where e is the energy density and n is the positive 
index, we have obtained exact solutions. The effect of the bulk viscosity is 
to produce a change in the perfect fluid. 
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